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The current fleet of large commercial aircraft has successfully achieved FAA noise certi- 
fications because of, in part, the successful application of uniform passive duct liner treat- 
ments to control engine system noise. One goal of NASA’s engine system noise reduction 
program is to develop technologies to improve the sound absorbing properties of duct liner 
treatments so that they remain effective in modern turbo fan engines. One such technology 
being studied is checkerboard or periodic axially and circumferentially segmented liners. 
A preliminary assessment of the potential of this technology was conducted by applying 
uncertainties associated with manufacturing, installation, source structure, and tonal fre- 
quency to a liner developed using deterministic design methods and generating a measure 
of improvement with respect to a uniform liner subjected to the same uncertainties. De- 
terministic design and analysis of the candidate checkerboard liner showed that it obtains 
a 1.5 dB per duct aspect ratio improvement in liner attenuation over a similarly designed 
uniform liner. When uncertainties in liner impedances, source structure, and frequency 
are considered, the performance of the checkerboard liner drops off dramatically. The final 
results of this paper show that the candidate checkerboard liner has a less than 25 percent 
chance of outperforming the uniform liner when moderate levels of uncertainty are con- 
sidered. It is important to note that this study did not include the effects of mean flow 
on liner performance and, more important to note, that as a gradient based optimization 
process was used to design the checkerboard liner, it is almost certain that a global optimal 
design was not found for the candidate checkerboard liner. Had it been possible to find a 
better deterministically performing checkerboard liner, the probability that this candidate 
liner would outperform the uniform liner would certainly have been higher. 


Nomenclature 

A nrn = magnitude of acoustic modes, Pa 

c = sound speed in duct, m/s 

F = liner attenuation, dB 

/ = frequency, Hz 

H,W,L = height, width, and length, m 

i = \/-T 

k = (cj/c), freespace wavenumber, m -1 

n = unit normal vector 

p = time independent acoustic pressure, Pa 
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probability 

normalized acoustic resistance and reactance 
dimensional time, s 

transverse, spanwise, and axial coordinates 
R + iXi normalized acoustic impedance 
three-dimensional Laplace operator, m -2 
increment of attenuation, dB 
phase of acoustic modes, radians 
ambient density, kg/m 3 
characteristic impedance of air, kg s -1 m -2 
27 r/, angular frequency, s -1 

exit plane 

uniform and checkerboard liner quantity 
the Jth impedance segment 
mode order in cross duct directions 
source plane 


I. Introduction 

The concept of the checkerboard liner is generally attributed to Mani, 1 who having shown experimentally 
that a circumferentially segmented liner could provide better attenuation than a uniform liner, proposed that 
a liner that is segmented both axially and circumferentially (i.e. , the checkerboard liner) may have potential 
for additional improvements. The paucity of research on the sound attenuation produced by checkerboard 
liners is due to a lack of fully 3-D analysis codes for modeling the liner attenuation within the duct. To 
date, no experimental and only four analytical studies 2-5 of sound attenuation produced by checkerboard 
liners have been reported. Howe 2 compared sound attenuation produced by uniform, axially segmented, 
circumferentially segmented, and checkerboard liners. The primary conclusion of that study was that the 
checkerboard liner offers the best means of increasing the attenuation over and above that of a uniformly lined 
duct. The fact that Howe 2 used a very restricted theoretical analysis (the kinetic theory of gases) to model 
liner attenuation, however, clouded the results of this study, and, Howe did not have liner optimization studies 
to support the conclusions of his paper. Watson et al. 3 conducted a study of three optimization techniques 
(iterative contour deformation, Davidon-Fletcher-Powell, genetic algorithm) for optimizing the attenuation 
produced by checkerboard liners. It was concluded that the genetic algorithm was the most reliable of the 
optimization techniques due to its ability to handle multiple local optima generated by the checkerboard 
liner. In a subsequent study, Watson et al. 4 combined the genetic algorithm with a numerical solution to the 
Helmholtz equation to demonstrate that an optimized 8-segment checkerboard liner gives considerably more 
attenuation than an optimum uniform or 4-segmented checkerboard liner. However, central processor unit 
(CPU) time and memory constraints and the occurrence of multiple local optima limited these results to a 
single frequency. Watson et al. 5 extended this analysis to multiple frequencies for the same configuration. 
The primary conclusion of the study is that an optimized checkerboard liner attenuates substantially more 
sound than an optimized uniform liner, and it has a broader off-design performance. Relative to the uniform 
liner, the checkerboard liner tends to be more effective at the higher frequencies. 

During the design optimization and assessment studies of Watson et al., 3-5 they observed that the 
checkerboard liner appeared significantly more sensitive to variations in the impedance values than did the 
uniform. This trend was also observed in the initial studies involving axially segmented, circumferentially 
segmented, and checkerboard liners. This is significant in that the values reported in many of the studies 
for liner impedances had significant digits greater than the manufactured tolerances of this type of liner. It 
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was also apparent that the structure of the source model strongly influenced the resulting distribution of 
impedances. This was most evident in the 2- and 4-segmented liner studies in which it was anticipated that 
some mirror image liner designs would appear. An investigation into their absence concluded that the source 
mode structure skewed the anticipated results. A further review of the literature found this observation had 
been previously reported. 

In the 70s, the axially segmented liner was studied extensively. Motsinger, 6 Kraft, 7 and Lester 8 performed 
some studies of particular interest to this investigation. However, the axially segmented lining configuration 
was not implemented in the current generation of aircraft engines because their optimum design requires 
an accurate description of the modal distribution of acoustic energy in the sound source. The problem of 
determining the modal amplitude and relative phasing of source modes within wide chord fan engine ducts 
is extremely difficult. 

Baumeister, 9 who studied the performance of multielement liners, provided a further indication of the 
dependence of optimum liner impedances on source structure. His observations were similar to those of 
Watson, concluding that multielement liners are most advantageous at high frequencies and did, theoretically, 
outperform uniform liners. They also observed that multielement liners suffered from large degradation in 
performance in response to changes in the assumed input source structure. 

The work presented in this paper was motivated by the need to provide an assessment of what performance 
gains may be attainable in an actual experiment or in service conditions by using checkerboard liners. The 
primary uncertainty that is the focus of this paper is the uncertainty in liner impedance associated with the 
manufacturing and installation of these liners in scale model engines and test facilities. Previous studies of 
the performance of checkerboard liners did not take into account these uncertainties. The level of variability 
in liner impedances is taken to be plus or minus one-tenth of a pc, about the nominal liner impedance value 
for this investigation. This value is more of a goal of NASA’s engine noise reduction systems program than a 
true reflection of the current state of the art. Previous publications 9 have suggested a variance in impedance 
of approximately twice this value. Other sources of uncertainty associated with operational conditions are 
assumed to be independent of liner uncertainties and to be taken into account through applying uncertainties 
to the source structure and frequency. The measure of the influence of these uncertainties on the performance 
of checkerboard liners is presented with respect to the performance of a uniform liner subjected to the same 
level of uncertainties in manufacturing and installation, operating conditions, and source structure. The 
remainder of this paper is broken into sections covering various aspects of this assessment. The next section 
of this paper discusses the problem formulation and solution methodology. The third section will discuss the 
uncertain variables, assessment metric, and procedure. The final section will present results with respect to 
specific variables. A brief discussion of the results and references will follow. Thus, this study will assess 
the advantages and disadvantages of a periodic checkerboard liner configuration over the currently installed 
uniform liners with uncertainties. 


II. Problem Formulation 

Both the uniform and checkerboard liners are evaluated using a geometry corresponding to that of the 
grazing incidence tube in the NASA Langley Research Center Flow Impedance Test Facility. This geometry 
and the selected frequency range were chosen so that experimental tests could be conducted at a later date 
to corroborate the computed results. Figure 1 gives a schematic of the 3-D, zero flow duct, and right-handed 
coordinate system used to predict the liner attenuation properties. The checkerboard liner is axially centered 
between the source and exit plane. The upper and two side walls of the duct are rigid. Acoustic waves are 
propagated from left to right across the surface of the liner and into a termination section of duct that is 
designed to be nonreflecting to plane acoustic waves. The test section is 0.812 m long and has a cross section 
of 0.051 m x 0.051 m. The source and exit planes are located in hard wall sections of the duct, 0.203 m from 
the leading and trailing edges of the liner, respectively. 
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A. Governing Equations, Solution, and Objective 


The mathematical problem is to find the 
solution to the Helmholtz equation sub- 
ject to the liner impedance boundary con- 
dition along the periphery of the duct, 
the exit impedance boundary condition 
at the duct exit, and a known sound 
source condition at the source plane (z = 
0). The Helmholtz equation is 

V 2 p + k 2 p = 0 (1) 


where a time convention of the form e luJt 
is assumed. The boundary condition 
along the rigid-wall portions of the duct 
is equivalent to the requirement that the 
gradient of acoustic pressure normal to 
the wall vanishes, that is 



( 2 ) 



pressure, p 6 


Figure 1. Schematic of grazing incidences impedance tube and 
liner configuration. 


where n is the unit normal vector to the duct wall and is taken as positive when pointing into the fluid. 
At the duct exit (z = L), the ratio of acoustic pressure to the axial component of acoustic particle velocity 
must equal the exit acoustic impedance. When expressed in terms of the acoustic pressure, this boundary 
condition is 

dp = -ikp 

dz Cexit 

A dimensional termination impedance of pc (£ e = 1) is used as a reasonable proxy impedance for all 
modes. It is not possible to know a priori the resultant impedance of some arbitrary superposition of modes 
incident upon the computational boundary plane. This exit impedance will not lead to an exact simulation 
of the physical termination of the grazing incidence tube because it may not allow higher order modes to 
transmit through the computational boundary properly. However, it is a more accurate representation of 
the grazing incidence tube termination than a nonreflecting boundary condition. 

The wall lining is assumed to be locally reacting, so that the liner boundary condition is 


dp ikp 
dn ( 


At the source plane (z = 0) of the duct, the sound source pressure, p s , is known. Therefore, the sound 
source boundary condition is 

P = Ps (5) 

Eqs. (1-5) form a well posed boundary value problem that one can solve to determine uniquely the 
acoustic pressure field in the duct. An exact analytical solution for this field is not available for a general 
input data set; therefore, a numerical method is used to obtain the solution for this field. 

In previous publications, this system of equations was solved using a research grade finite element pro- 
gram. This program utilized the scientific library on Silicon Graphics, Inc. computers. In this investigation 
the solution methodology was transferred to the commercial code Actran. 10 This code directly calculated 
the radiated power out the end of the finite length duct. A duct with rigid walls was used as a reference 
point from which to calculate the objective function. 
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Watson et al. 4 used the reduction in sound power from the source to the exit of the duct as their objective 
function. This function was maximized to determine the optimal attenuation of the liner. This investigation 
also uses this measure of the reduction in sound power as the primary measure of performance for individual 
liners. This function is normalized as before with respect to the duct and liner dimensions and has units of 
decibels per (L 2 —Li)/H. The radiated sound power is calculated directly by the analysis program Actran 10 
and no further post processing is necessary other than normalizing the result. 

III. Uncertainty assessment 

The desired metric for the assessment of the liner performance under uncertainty can be stated as the 
probability that the attenuation of the checkerboard liner is greater than the attenuation of the uniform 
liner by an arbitrary increment A. This increment, along with the probability of attaining it, can be used to 
assess the merits of further exploring this concept subjectively. Mathematically this probability is stated as 

p (F c ((i,Q nm ,u) > F u (( u ,O nm ,Uj) + A) (6) 

where , ( u , 0 nm , and uj are independent random variables quantifying uncertainties in the impedance of the 
I th liner segment, impedance of the uniform liner, modal source phase, and source frequency, respectively. 
This paper discusses each of these separately in successive sections. The probability stated in Eq. (6) is 
calculated from the statistics of the uniform and checkerboard liner attenuations as 

/ +oo r+oo 

/ pdf(F c ,F u ) dF c dF u (7) 

- oo J F u A 

where pdf(F c , F u ) is the joint probability density of the liner attenuations. Since F c and F u are independent, 
Eq. (7) can be written as 

/ +oo 

[1 - cdf(F u + A)] pdf(F u ) dF u (8) 

-OO 

where cdf(F u + A) is the cumulative distribution function for the checkerboard liner evaluated at F c = 
F u + A. The probability and cumulative probability distribution functions are obtained by direct Monte 
Carlo simulation, and the integration in Eq. (8) is carried out numerically. 

A. Modal Source Model with Uncertainties 

A modal pressure distribution appropriate for a multimodal noise source located in a rigid-wall section of 
duct, 

N M 

/'(•'•• y, <>) = EE A n m^ dnrn cos (mry/W) cos (rmrx/H) (9) 

n = 0 m = 0 

is used to model the sound source. The upper limits of summation are only taken over the cut-on modes, 
which were assumed to have unit amplitude and zero phase (i.e., A nrn = 1 and 0 nrn = 0) in previous 
investigations. 5 * In this investigation, the amplitude, A nm , is again assumed to have unit amplitude for all 
modes; however, more of the phase of each mode, # nm , is considered uncertain. The phase will be allowed 
to vary from —it to i r with all angles being equallikely. Note the use of Eq. (9) in this and the previous 
investigations contains an assumption that the rigid-wall duct source does not vary appreciably with the 
introduction of the liner into the duct. 
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B. Liner Impedances with Uncertainties 


The uncertainty in liner impedance arises primarily from manufacturing and installation considerations. 
Previous investigations that have considered uncertainty have used values of 0.2 pc. This investigation uses 
the value of 0.1 pc. This value is not representative of the current state of manufacturability but is a reflection 
of the goal of NASA’s engine noise reduction systems program. As in previous publications, the resistance 
and reactance of the liner segments will be assumed to be independent random variables. The mean value 
for these variables will be the values that were obtained from the deterministic optimization study performed 
at 7 kHz. These results will be revisited and presented in the first section of the result section of this paper. 

C. Frequency Dependence and Uncertainties 

The variation in frequency arises quite naturally in experimental and real world situations. The magnitude 
of variation of frequency is assumed only to be a few percent. In this investigation, the mean frequency is 
that for which the liner impedances were optimized. Thus, the the frequency is assumed to be Gaussian 
distributed with the mean value of 7 kHz and a standard deviation of 180 Hz. This frequency range was also 
chosen such that the structure of the source would not change over this range. That is, no modes would be 
cut-on or cut off in this frequency range. The frequency dependence of the liner impedances must also be 
modeled. Watson et al. 5 used a suitable model for this facility in their previous investigation. This model 
takes the form 


= Rj + i\i = Ri T i cot (kdj) (10) 

where Ri and \i are the resistance and the reactance, respectively and are taken to be specified at the 
design frequency. The resistance is assumed independent of frequency, and the reactance is assumed to vary 
as the cotangent of frequency. The value dj in Eq. (10) is the depth of the liner. This model is consistent 
with the experimental findings of Gallman and Kunze 11 for a similar facility. 

D. Monte Carlo Technique 

A Monte Carlo technique was used to estimate the probability and cumulative probability density functions 
of the liner attenuation. The commercial software package iSIGHT 12 was used for this process. With 
variations in impedances, frequency, and modal phase, the distribution of the uniform liner attenuation was 
estimated using approximately 400 samples. The estimate for the segmented liner with four independent 
impedances, as well as the uncertainties in frequency and phase, used 1000 samples. The criteria for the 
number of samples was subjective in that only the minimum number of samples to accurately capture the 
distribution functions for each of the individual variables were generated. The convergence of the mean and 
variance of the attenuations were, however, verified. This investigation did not apply formal methods to the 
sampling. The distribution functions were estimated from the Monte Carlo results using standard histogram 
methods with 50 equally spaced bins. These results will also be presented in the second section of the results 
section of this paper. 


IV. Results 

The result section of this paper is broken into three sections. The first section will discuss the deterministic 
optimization results and examine briefly the characteristics of the acoustic field inside the lined ducts with 
respect to the unlined duct. The second section will discuss the results from the Monte Carlo simulation using 
the liner impedances, source phase, and frequency as sources of uncertainty. The final section will present 
the results from the assessment and present the trade-off between liner performance and the likelihood of 
obtaining the performance. For clarity the geometry of the duct and liner arrangement is repeated here. 
Figure 1 gives a schematic of the 3-D, zero flow duct, and right-handed coordinate system. The checkerboard 
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liner is axially centered between the source and exit plane. The acoustic waves are propagated from left to 
right across the surface of the liner. The test section is 0.812 m long and has a cross section of 0.051 m x 
0.051 m. The source and exit planes are located in hard wall sections of duct, 0.203 m from the leading and 
trailing edges of the liner, respectively. The radiated sound power is calculated at the exit plane. The cut-on 
frequencies for the rigid-wall duct modes of interest and the next higher modes are presented in Table 1 


Table 1. Mode Cut-On Frequencies 


Mode Number (n,m) 

Frequency (hertz) 

(0,0) 

N/A 

(1,0) 

3346.457 

(0,1) 

3346.457 

(1,1) 

4732.602 

(0,2) 

6692.914 

(2,0) 

6692.914 

(1,2) 

7482.905 

(2,1) 

7482.905 


The orientation of the impedances 
in the checkerboard liner are shown in 
Fig. 2. Impedances (i and £2 are clos- 
est to the source plane. The numerical 
model of the duct is comprised of 90,457 
nodes using 79,872 linear hexahedral fi- 
nite elements. There are 17 finite element 
nodes in the cross duct directions and 313 
nodes along the length of the duct. This 
is the identical numerical model used in 
previous investigations. 

A. Optimization Results 

This paper presents limited optimization 
results with respect to those presented in 
previous publications. During the prepa- 
ration of this paper and conduction of the 
studies outlined, an error in the source 
structure of the previous publications was 
detected. There were insufficient time 
and resources available to rerun all of the 
deterministic optimization results using a 
genetic algorithm. The deterministic op- 
timization procedure used in this inves- 
tigation was gradient based and, subse- 
quently, could not be guaranteed to find 
global optima particularly in the case of 
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Figure 2. Schematic of uniform and periodic checkerboard liners. 


the checkerboard liner. The impedances tabulated in the previous studies were, however, used as starting 
points for the optimization in this study. Table 2 presents the impedance values obtained in this investigation 
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at design frequency of 7 kHz with all modes of equal amplitude and zero relative phase. In this table the 
uniform liner impedance is ( u and C15C25C3 an d (4 are the impedances of the checkerboard liner. 


Table 2. Impedance Values Prom Optimization 


Segment 

impedance pc 

liner depth(mm) 

Cu 

2.37 — 2.35i 

3.09 

Ci 

3.46 - 3.7(k 

2.04 

C 2 

0.62 - 3.25 i 

2.31 

Cs 

2.25 - 0.42z 

9.04 

C4 

1.48-0.8B 

6.88 


The performances obtained by this process were 2.28 dB per duct aspect ratio for the uniform liner and 
3.74 dB per duct aspect ratio for the checkerboard liner. Table 2 also presents values for the liner depth 
obtained from Eq. (10). No attempt was made to constrain these values to physically realizable quantities. 
It is clear from this table that most of these liners are nonphysical, having depths less than the thickness 
of typical liner face sheets. This practice, however, is consistent with previous related publications. The 
process of applying such a constraint would not have significantly complicated the design and assessment 
process. The different acoustic fields inside the uniformly, checkerboard, and unlined ducts are shown in 
Fig. 3, respectively. The unlined duct is presented for completeness and does show that the acoustic field is 
correctly resolved. It is also important to note the strength of reflections off the termination. 



(a) No liner. 


(b) Uniform liner. 


(c) Checkerboard liner. 


Figure 3. Comparison of real part of pressure in duct with no liner, uniform liner, and checkerboard liner at 
7 kHz. 


B. Monte Carlo Results 

Monte Carlo simulations were conducted on the uniformly and checkerboard lined ducts about the deter- 
ministic optimum solution. The uncertainties used were those previously stated for liner impedance, source 
phase, and frequency. For the uniform duct, this resulted in eight uncertain variables; and for the checker- 
board lined duct, 14 uncertain variables. This investigation used 1000 random samples for the checkerboard 
liner and 423 random samples for the uniformly lined duct. The odd number of samples for the uniformly 
lined duct is a product of a premature termination to the run. The simulations did appear to converge; 
however, no formal methods were applied to the sampling process. For completeness the distributions of the 
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uncertain variables are shown in Fig. 4 for the uniformly lined duct and in Fig. 5 for select variables of the 
checkerboard lined duct. It is important to note that the frequency, being Gaussian distributed, does span 
cutoff frequencies of several modes both above and below the design frequency. 



(a) Resistance Distribution. 


(b) Reactance Distribution. 




(c) Frequency Distribution. (d) First Mode Phase Distribution. 


Figure 4. Probability distribution functions of variables for uniform liner Monte Carlo simulation. 


The resulting probability density functions for the uniform and checkerboard liners are shown in Fig. 6. 
Note that the attenuations are negative in this figure. Larger negative numbers indicate greater attenuation. 
This is the opposite of the deterministic results. The figures show that the variance in the uniformly lined duct 
is smaller than that for the checkerboard lined duct. This is consistent with previous findings and indicates 
that the uniformly lined duct is less sensitive to uncertainties than is the checkerboard lined duct. The 
distribution for the checkerboard lined duct is bounded on the upper limit of attenuation by the deterministic 
value of —3.74 where as the distribution of the uniformly lined duct is not. The deterministic value, —2.28, 
of the uniformly lined duct falls nearer the mean of the distribution. This result was unexpected and had 
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(a) First Segment Resistance Distribution. 


(b) First Segment Reactance Distribution. 



Frequency 



(c) Frequency Distribution. 


(d) First Mode Phase Distribution. 


Figure 5. Probability distribution functions of select variables for segmented liner Monte Carlo simulation. 
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not been seen until uncertainties in frequency and phase were included in the simulations. An examination 
of the raw data yielded little information or insight into the cause of the resulting distribution, however, 
the resulting distribution appears to be comprised of at least two distinct effects combining independently. 
The first appears to be a distribution comprising the right lobe of the resulting distribution and believed 
primarily due to uncertainties in source structure. The second is believed to be responsible for the left 
lobe and arising from the frequency distribution. The frequency distribution spans a range in which duct 
modes are being cut-on and cutoff and also in which the frequency versus liner attenuation curve may have 
a steep slope. This second effect is believed responsible for the attenuations greater than the deterministic 
attenuation. This discussion of the distribution is speculative but is consistent with anecdotal information. 

The distribution of the attenuation of the checkerboard lined duct is characteristic of a highly tuned dy- 
namic system. Performance drops off dramatically with even the slightest change in the operating conditions. 
This figure indicates that there is virtually no possibility of obtaining the 3.74 dB per L/H of additional 
attenuation predicted by the deterministic analysis. This observation was consistent with previous studies 
and has been one reason why such segmented liners have not been considered viable. 
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(a) Uniform liner. 


(b) Checkerboard liner. 


Figure 6. Probability distribution functions of duct attenuation for checkerboard and uniform liners. 


C. Liner Performance with Uncertainty 

The probability that the checkerboard liner would outperform the uniform liner by an increment A is plotted 
in Fig. 7. The first and obvious thing to note is that this checkerboard liner has a less than 25 percent chance 
of outperforming the uniform liner under the range of conditions specified in the Monte Carlo simulations. 
The deterministic performance gain of approximately 1.5 dB per duct aspect ratio is virtually nonexistent 
when the influence of uncertainties are included in the assessment. This result is valid only for this liner 
configuration with these values of impedance at this frequency. The results from a deterministically better 
designed checkerboard liner could be substantially better. This result does not address the full performance 
of the liner over the full frequency range of previous investigations and, therefore, does not address whether 
the uniform liner or the checkerboard liner has better broadband characteristics. 

V. Discussion 
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The results of 
this investigation show 
the potential ad- 
vantages and dis- 
advantages of a checker- 
board liner. The 
configuration of the 
checkerboard liner 
was determined by 
deterministic opti- 
mization and the 
assessment of its 
performance carried 
out by Monte Carlo 
simulation. The 
liner design vari- 
ables were those of 
the DC component 
of flow resistance 
and liner depth. 

This investigation 
did not constrain 
these quantities to 
be physically real- 
izable. This inves- 
tigation, also, did 
not take into ac- 
count the influence 

of the flow field and other characteristic dimensions of locally reacting duct liners. Results of this investiga- 
tion did show that this checkerboard liner was unlikely to outperform a uniform liner under the assumptions 
of uncertainty in impedance, source phase, and frequency oscillations. 
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Figure 7. Probability of obtaining more attenuation from checkerboard liner over uni- 
form liner. 
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